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Abstract. Quantum cellular automata (QCA) are models of quantum computation of 
particular interest from the point of view of quantum simulation. Quantum lattice gas 
automata (QLGA - equivalently partitioned quantum cellular automata) represent an inter- 
esting subclass of QCA. Prior work on QCA has investigated the relationship between these 
two classes of models. In the present paper we establish necessary and sufficient conditions 
for unbounded, finite Quantum Cellular Automata (QCA) (finitely many active cells in a 
quiescent background) to be Quantum Lattice Gas Automata. We define a local condi- 
<**j . tion that classifies those QCA that are QLGA, and we show that there are QCA that are 

not QLGA. We use a number of tools from functional analysis of separable Hilbert spaces 
and representation theory of associative algebras that enable us to treat QCA on finite but 
unbounded configurations in full detail. 



1. Introduction 

Feynman first noted that simulating the full time evolution of quantum systems on classical 
computers is a hard problem, and that one might use one quantum system to efficiently 
emulate another [14-16]. Feynman's suggestion has been subsequently developed in the 
field of quantum simulation, which has attracted considerable theoretical and experimental 
attention in recent years [5,27,30,31,40,41,48,54,56,58]. 

Simulation of quantum systems on current classical computer hardware is a well estab- 
lished field with simulations of diverse systems from quantum chemistry to the structure 
of the proton. For large systems these simulations rely on approximate methods, such as 
semiclassical treatments or the Monte Carlo method, which scale only polynomially in the 
problem size. The results of even approximate methods are of interest as they address issues 
which are not accessible in any other way. Classical simulation of quantum systems will 
remain a hard problem for decades to come, and one may expect useful quantum computers 
to appear on this timescale. 

Amongst both classical and quantum simulation methods, cellular automata, lattice gas 
and random walk methods can be singled out for their simplicity. In 1948 von Neumann 
set out to show that complex phenomena can arise out of many simple, identical interacting 
entities. Following a suggestion by Ulam he adopted an approach in which space, time 
and the dynamical variables are all discrete. The result was the cellular automaton, a 
homogeneous array of cells with a finite number of states evolving in discrete time according 
to a uniform local transition rule [52]. 
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Figure 1. Non-trivial collisions of the HPP lattice gas: The dynamics takes 
place on a cartesian lattice in two dimensions. A pair of incoming particles 
with opposite velocities in the North-South direction become a pair of particles 
with opposite velocities in the East- West direction, and vice versa. 

In the context of physical simulation, classical cellular automata led to lattice-gas models 
- also called partitioned cellular automata - which can simulate diffusive processes and fluid 
mechanics for both simple and complex fluids [44] . Lattice gases are the only broad class of 
cellular automata models that have enjoyed wide success in quantitative modeling of physical 
phenomena. 

In the quantum setting, quantum cellular automata [2,3,13,33,45,50,55], quantum lat- 
tice gases [8-10,34-38] and quantum random walks [1,12,28] have all attracted considerable 
attention. In particular the one dimensional cellular automata for the Dirac equation, orig- 
inally described by Feynman in a problem in [17], has been investigated and generalized by 
Meyer into the quantum lattice gas model [34-38]. Similar models were investigated inde- 
pendently by Boghosian and Taylor [8-10]. Recently, optical implementations of quantum 
random walks have demonstrated topologically protected bound states, [29] and simulations 
of the Dirac equation have been performed in trapped ions [20]. 

1.1. Classical lattice-gas cellular automata. Perhaps the most widely explored appli- 
cation of cellular automata has been in the field of fluid dynamics. They have been used 
extensively for modeling hydrodynamics since Frisch, Hasslacher, and Pomeau [19], and Wol- 
fram [57] showed that it is possible to simulate the incompressible Navier-Stokes equations 
using discrete boolean elements on a lattice. The update rule for all lattice gases takes place 
in two steps, propagation and collision. For the simplest such models one represents the 
presence or absence of a particle by a bit. Each bit at a lattice site is associated with a link 
in a lattice and a corresponding velocity. During propagation the bits move in the direction 
of their velocity vectors along the links, retaining their velocities as they do so. During 
the collision step, the newly propagated bits are modified according to a purely local rule. 
Specification of collision rules completes the description of a lattice gas automata (LGA). If 
the collision process is stochastic, a single particle in such a model will follow a random walk, 
and so one may regard classical lattice gas models as equivalent to classical multiparticle 
random walks with an exclusion principle. 
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The simplest two-dimensional example of a lattice gas model is the HPP model [24] . Here, 
the underlying lattice is cartesian, hence there are four lattice vectors per site, corresponding 
to a particle moving north, east, south, or west. There may be at most one particle per 
direction per site, leading to at most four particles at each lattice site. The non-trivial 
collision of the HPP gas is shown in Figure [TJ 

The distinction between cellular automata and lattice gases has been the source of some 
debate [25,49]. A homogeneous cellular automata utilizes the same neighborhood and the 
same update rule at every site and every timestep, as in the Game of Life. The lattice- 
gas evolves in two distinct substeps: propagation and collision (described in detail below). 
The connection between the models may be made by introducing the concept of a block- or 
partitioned- cellular automata, which alternates between two update rules and two neigh- 
borhoods. The simplest of these automata employ a Margolus neighborhood [32]. As an 
example, we shall construct the two rules and two neighborhoods which, when alternately 
applied to a cellular automata, yield the HPP lattice gas dynamics described above. The 
update rules are shown in Figure [2J The update rules act on a block of four cells, which we 
shall refer to as upper left (UL), upper right (UR), lower left (LL), lower right (LR). The 
propagation rule is applied identically to all states: the occupation of each cell in the four 
cell block is exchanged with the occupation of the cell diagonally opposite. The collision 
step is identical to the propagation step except when UR and LL are the only cells occupied 
or when LR and UL are the only cells occupied, in which case the occupancies of these pairs 
of cells are exchanged. Graphically this corresponds to either diagonal of the four cells being 
occupied. In these cases the state is flipped to the other diagonal occupation, as shown in 
Figure [21 
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Figure 2. Update rules for the HPP partitioned cellular automata. The top 
two boxes show the collision rules in which the configurations of the cells which 
are modified in a non-trivial way. The bottom boxes show the update rule 
which implements propagation in the automata. The bottom rule is applied 
to both neighborhoods of the automata. A, B, C, D represent any occupation 
state of the four cells, and the propagation update is the permutation of these 
occupations shown. 

Clearly, if either of these rules was applied to a constant neighborhood the resulting 
dynamics would be trivial. Neither rule can change the occupation outside the four cell 
block and both rules are self inverse, so applying them twice to the same block returns the 
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original configuration. In our partitioned automata we specify a different neighborhood for 
each rule in order to obtain the HPP dynamics. The Margolus neighborhoods, the cells of 
the automata and their relationship to the original cartesian lattice on which we defined 
HPP model are shown in Figure |3J The two neighborhoods are the two ways of partitioning 
a cartesian grid into square tiles of four cells. The original cartesian lattice on which the 
HPP lattice gas was defined is the lattice of diagonals with a lattice node at the center of 
every collision neighborhood. The lattice gas vertices are shown by open circles in the third 
Figure in Figure [3j 



Figure 3. Neighborhoods for collision (left) and propagation (middle). The 
boundaries of the four cell blocks on which the rules shown in Figure |2] are ap- 
plied are shown by the bold lines. The figure on the right shows the boundaries 
for the collision neighborhood. The central figure shows the boundaries for the 
propagation neighborhood. The figure on the right shows the correspondence 
between the cells of the partitioned automata and the lattice of the lattice gas. 
Bold lines denote the directions along which lattice gas particles move. Open 
circles denote vertices of the lattice-gas lattice. 

In this approach, the original cells of the automata are grouped together into blocks of four 
to make the sites of the HPP lattice-gas model. The partitioned dynamics is reinterpreted 
- one step of the dynamics (which acts on one partitioning of the original cells) now acts 
locally on a single lattice-gas site. This local action is termed the collision or scattering step. 
The second step of the dynamics (which acts on the other partitioning of the original cells 
in the Margolus scheme [32]) implements propagation of information between the lattice-gas 
sites. In the case of the HPP model, where the alphabet of the original cells is {0, 1}, this 
step has the interpretation of particles moving from one lattice gas site to another. The four 
sub-cells of a lattice gas site are then reinterpreted as lattice vectors, and in the propagation 
step lattice-gas particles move in the direction of their lattice vector to a new site. Hence we 
have taken a partitioned cellular automata and reverse-engineered a lattice-gas model from 
it. 

The general question of when a classical cellular automata can be reinterpreted in this way 
was recently taken up in [49]. Because we know, from the counter example of [3], that there 
are QCA that are not partitioned CA, we shall proceed by first defining classical lattice-gas 
automata and proving a general condition for when a QCA is a QLGA. In the process we 
shall extend our definition of a classical lattice gas to the quantum case. 

Consider the lattice to be Z n . Each point x on this lattice is a site of the lattice gas and 
has the same neighborhood Af x . We may consider each neighborhood M x to be a translation 
(by x) of the neighborhood M of the site x = 0, where M is a finite subset of Z n . Hence 
the neighborhood sites of site x are simply x + z z E Af. The state of a single site may 
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be constructed by assigning a substate to each of a set of lattice vectors - each vector 
corresponding to a neighborhood element z G M . The sub-state k x (z) of each lattice vector 
z at lattice site x takes values in an alphabet C z . For the simplest lattice gases, each lattice 
vector may have either one or zero particles present and so C z = {0, 1} V z. However, more 
sophisticated models allow multiple types of particles, and three-dimensional models allow 
some vectors to carry more than one particle [7] , hence in general C z may contain arbitrarily 
many symbols and may vary with z. The state k x of a lattice-gas site x is therefore an 
element of the Cartesian product of all of the C z : 

k x = l[k x (z) e B = ]JC Z (1) 

The state of the lattice gas is an element of the set of infinite sequences S over B. 

We may now define the two substeps of the lattice gas dynamics. Firstly, the sub-states 
propagate to the lattice sites in the neighborhood: 

a : S 

n Hk x (z)^ n ]jk +z (z) 

Secondly, in the collision step, the state of each lattice-gas site changes locally under a map 
0: 

e : n c * — ► n c *- ( 2 ) 

With these classical definitions in hand, we now turn to the definition of quantum cellular 
automata. 

1.2. Watrous Quantum Cellular Automata. QCA have been studied, in considerable 
depth, by several authors. The main idea is to develop quantum models that retain the 
key features of a classical cellular automata: translation invariance and discreteness in time, 
space and dynamical variables. A classical deterministic cellular automata is defined as 
follows: 

(i) A lattice L, whose sites we shall label by x. 

(ii) A set of cell values Q with a distinguished quiescent state q . 

(iii) A neighborhood M . This translates to a neighborhood M x for each cell x, consisting 
of a set of cells as follows: M x = {x + z\z G A/"}. 

(iv) A configuration of the automata is a map a : L — > Q, If L is infinite then only a 
finite number of cells take values in Q that are not quiescent. 

(v) A local transition rule S which updates the value s G Q at x according to the values 
Si G Q of the cells in the neighborhood 5 : >->■ Q, subject to the condition that 

5 : gf 1 h. g . (3) 

Early work produced an extension of the classical cellular automata (CA) defined above 
based on local rules [2,13,33,50]. To avoid confusion with the definition of a Quantum 
Cellular Automata (QCA) that we shall introduce later we follow [45] and call these Watrous 
Quantum Cellular Automata (WQCA). A WQCA is defined as follows: 

(i) A lattice L, whose sites we shall label by x. 

(ii) A set of cell values Q with a distinguished quiescent state q . 
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(iii) A neighborhood M. This translates to a neighborhood M x for each cell x, consisting 
of a set of cells as follows: M x = {x + z\z G Af}. 

(iv) A state of the automaton, is a superposition of classical configurations. 

(v) A local evolution rule / which maps an element of to an element of W = C[Q], 

f:Q m ^ C[Q] 

subject to the constraint that if the neighborhood state is q^ then / assigns ampli- 
tude 1 to the quiescent state and amplitude to all other states. 

From the above data, one can construct a global evolution rule R by jointly evolving the 
state of each cell % in any basis state by the local rule / and extending by linearity. The 
global evolution R is required to be: 

(i) Unitary: R)R = I. 

(ii) Translation invariant: Let t z , for some z G Z, be the translation operator defined 
by its action on a classical configuration c = Yljez c j f° r c j Q V T- 

t z ■ II c i ^ II c i+ z 

t z is extended by linearity. We denote the group of translations: T = {r z } ze z- A 
linear operator R is translation-invariant if t z Rt~ 1 = R for all zGZ. 

By definition, the above model of WQCA is translation-invariant provided the rule / is 
the same for all cells. Unitarity for this model has been extensively investigated [13,33,50]. 
This model of QCA was found [3], in some instances of the finite, unbounded QCA, to allow 
super-luminal signaling. Super-luminal signaling refers to the faster than light signaling 
which occurs when the state of a cell (restricted or reduced state, defined in the next section) 
can be affected by that of another cell which is an unbounded distance away from it. This 
is observed when certain classical CAs with finite neighborhood schemes are quantized (the 
local rule for classical CA is taken as that for the QCA). 

1.3. Axiomatic QCA. An axiomatic approach has been developed to overcome the prob- 
lem of super-luminal signaling, a non-local phenomenon, seen in the WQCA models [3,4,45, 
46]. The requirement of causality is imposed in the definition of a QCA, with roots in the 
ideas of [6,46], to preclude this non-locality [3,4,11,45]. 

The axiomatic model, first developed by Schumacher and Werner [45], builds in the re- 
quirement of causality from the start. Causality can be defined either in terms of states or 
in terms of operators. For states, causality is the requirement that the reduced state of a 
cell depends only upon the reduced state of its neighborhood at the previous timestep (with 
the caveat that quantizing a classical cellular automata may change the neighborhood, as 
discussed in Section I4.3p . For operators, causality states that an operator localized on a 
single cell at time t will be localized on the neighborhood of that cell at time t + 1. 

In this approach a QCA is defined to be a unitary, translation-invariant, *-homomorphism 
of a C*-algebra which satisfies a quasi-local condition by which the global homomorphism R 
is uniquely determined by a local homomorphism of tensor factor (one-site) algebras to the 
neighborhood. The main claim of [45] states that the structure of such QCA's can be given 
by ^isomorphisms obeying a quasi-local condition as follows. Once the cells are grouped into 
super-cells, such super-cells partition the lattice. Then there exists a shifted partition with 
super-cells that overlap with the original cells in a Margolus partitioning scheme [32]. The 
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claim is that there is a set of unitary operators U, V, each acting on a cell of the respective 
partition, can be applied in succession, to obtain the global evolution. Hence the claim of [45] 
is that imposing the requirement of causality on quantum cellular automata requires them 
to be quantum lattice gas automata. 

However, a counter example to the analysis of [45] was found by Arrighi et. al. in [3], and 
the axiomatic approach was revised in that paper. The axiomatic one-dimensional develop- 
ment of Arrighi et al [3] is more restricted and defined on a Hilbert space of unbounded, finite 
configurations (finitely many active cells in a quiescent background). This set of sequences, 
called the set of finite configurations, is embedded into some abstract Hilbert space in which 
the elements of this set comprise an orthonormal basis. This is a countable-dimensional 
Hilbert space - the Hilbert space of finite configurations. In this context, they define a 
QCA as a unitary, translation-invariant and causal operator on the Hilbert space of finite 
configurations. 

In the present work we further develop the axiomatic approach and place it on a new 
mathematical foundation. We incorporate in the development of axiomatic QCAs the topol- 
ogy of the Hilbert space. This is an essential part of any study of the algebra of operators, 
particularly infinite dimensional algebras. This is intimately connected, in this case, with the 
idea of local algebras and the part they play in the determination of the global evolution. 
This topology provides the foundation for connecting the local algebras to the algebra of 
bounded linear operators of the Hilbert space. With this topology, we can investigate the 
space of bounded linear operators as a von Neumann Algebra (besides being a C* algebra), 
making the local subalgebras a useful tool. 

We develop a model on the same space of finite but unbounded configurations as Arrighi 
et al [3]. To make a Hilbert space from it, we give an inner product structure to the set of 
finite configurations that is a natural extension of the inner product of the Hilbert space of 
each cell. We define the Hilbert space of finite configurations as the £ 2 completion of the 
linear span of the set of finite configurations. This defines the separable topology of the 
Hilbert space and the induced inner product is used to define the norm, weak, and strong 
and topologies of the bounded linear operators on the space. At the outset, we show how 
the subalgebra of local finite dimensional operators is connected to the algebra of bounded 
linear operators on the whole space (Theorem 13.71 the Density Theorem). This provides the 
background for the development of some of the key ideas, and underpins the proof of the 
Structural Reversibility Theorem (Arrighi et al [3]). 

The main result of this paper is that QLGA can be characterized by a local condition on 
QCA. This condition pertains to the set of image algebras under the global evolution of the 
neighboring cells. The condition for a QCA to be a QLGA is that the local pieces of these 
image algebras generate the full cell algebra. It is clear from our development that a central 
result, Theorem I3.10[ about the tensor product decomposition of Hilbert space of a single 
cell is not true in general. In particular, to have that decomposition, the aforementioned 
subalgebra condition is needed. The counter-example in two dimensions to Schumacher and 
Werner's [45] result, that Arrighi et al [3] have found, occurs precisely because that QCA 
violates this condition. This means that there is a class of QCAs that are not QLGAs and 
that require further exploration. 

Our proof methods differ from those in prior works significantly, in that we use the theory 
of infinite-dimensional von Neumann algebras, representations of associative algebras in finite 
dimensions, and the theory of separable Hilbert spaces and bounded linear operators. 
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The current paper is structured as follows. In Section [2] we introduce the Hilbert space of 
interest for our automata. In Section [3] we introduce the axiomatic definition of QCA and 
prove several properties of operators that are translation invariant or causal, or both. We 
extend our definitions of the propagation and collision operators to the quantum case and 
formally define a QLGA before proving our main theorem. In Section H] we give a number of 
examples of QCA, for the case when our condition is satisfied and when it is not. We close 
the paper with some conclusions and discussion. In the interests of making the present paper 
self contained we collect a number of background and corollary results in the appendices. 

2. The Hilbert space of finite configurations 

We begin with a finite set of symbols Q containing a special quiescent symbol qo, and an 
infinite lattice C = Z™ in n dimensions. Let us denote by W the Hilbert space of formal 
linear combinations of symbols in Q, i.e. W = C[Q] (in this paper all of the vector spaces 
will be over C). This Hilbert space W is the Hilbert space of quantum states of a single 
cell (for the rest of the paper, we will use the terms cell and site to mean an element of the 
lattice £). Let the dimension of W (cardinality of Q) be dimjy = dw The set of basis 
elements of W corresponding to the symbols in Q is denoted Bq\ 

B Q = {\q):qeQ} (4) 

W has an inner product: 

(a\0) = Y,*p9 (5) 

where a = £ 9eQ a 9 |g), = EgeQ^I?)- 

Informally, the Hilbert space of the automata is the infinite tensor product of copies of 
W. The idea of infinite tensor products was first studied by von Neumann in his landmark 
paper "On Infinite Direct Products" [51]. Later, Guichardet [23] introduced a more modern 
way of describing an infinite tensor product of Hilbert spaces enumerated by a countable set. 
This is sometimes referred to as the incomplete tensor product construction. From countably 
many copies of the same Hilbert space W the incomplete tensor product is obtained as the 
inductive limit of an ascending chain of finite tensor products of W. This space has a natural 
basis that we call the set of finite configurations. The Hilbert space of finite configurations is 
the incomplete tensor product construction of Guichardet with the set of finite configurations 
as its basis [23]. 

Definition 2.1. The set of finite configurations, denoted by C, is the set of simple tensor 
products {<g) x&n \c x ) : \c x ) G B Q , all but finite \c x ) = \q Q )}. 

Thus, C is a countably infinite set. Let |c) = (£) xeZ n \c x ), \d) = (£) xeZ n \c x ') G C. Define 
the inner product of the elements |c), |c'): 

(c\c) = Tl xeZ n (c x \c x '} 

and extend it by linearity to get an inner product on span(C) (the set of finite linear combi- 
nations of elements of C). 

Definition 2.2. The Hilbert space of finite configurations, denoted by He, is the £ 2 com- 
pletion of span(C) under the norm induced by the above inner product. C constitutes an 
orthonormal basis of He- 
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By definition, He is a separable Hilbert space since it has a countable orthonormal basis 
C (this follows from a standard theorem on Hilbert spaces in texts on analysis, for example, 
Folland [18], Proposition 5.29, pg. 176). 

Definition 2.3. The neighborhood M C Z n is some finite set. Then the neighborhood of a 
cell x <E Z n , denoted by Af x , is a translation of Af to x. 



The state of the QCA, p, is described by a density operator. 
Definition 2.4. A density operator, p, is a positive trace class operator on He, with tr(p) = 1. 



In this section we introduce the axiomatic definition of QCA in which the requirement 
of causality is built in from the start. We show that a QCA on the Hilbert space of finite 
configurations must possess a quiescent state as an eigenstate of its global evolution. We 
define local operators, and connect the algebraic structure of these operators to the topo- 
logical structure of operators on the Hilbert space of finite configurations. We then consider 
the action of the global evolution operator on the local operators and define the requirement 
of causality in terms of the evolution of local operators. This definition of causality in the 
Heisenberg picture was termed dual causality or structural reversibility by [45] and [3] . Struc- 
tural reversibility is useful because it means that the analysis of the causal global evolution 
can proceed by analysis of finite dimensional algebras of local operators. 

We continue from the definitions of Section [2J We need a few more definitions to properly 
state the requirements for an axiomatic model of QCA. For a finite subset D G Z™, we 
introduce the idea of a co-D space, which will be of use for a number of results. 

Definition 3.1. Let D G Z n be a finite subset. Define the set of co-.D configurations to 
be Cjj = {<S> ie z n \D \ c i) '■ c i e Qi au finite |cj) = l^o)}- Let the inner product on Cjy 
be induced by the inner product on W (0, as for C. Then co-D space, denoted : He—, is 
defined as the completion of span(C;p) under the induced £ 2 norm. 

Any operator E G B{Hc) can be expressed as a finite sum: 



where {e{\ C End((^) fc6Z) W) is a linearly independent set, and {E[} C B(Hc^)- 
We define the restricted or reduced density operator. 

Definition 3.2. Let p be a state (density operator) on He- The restriction of p to D is 
defined to be the partial trace of p over all the tensor factors corresponding to indices not 
in D: 



The partial trace can be defined as follows. We let {(u^)} be some orthonormal basis of 
the co-D space He— Then p can be expressed as: 



A4 = {k + x : k G Af} 



3. Axiomatic Definition of QCA 




(6) 



P\d = to{ x <£D}(p)- 



P = ^Px,y ® \w x )(w y 
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where p x>y G End((££) - eD Wj). The partial trace is defined as a sum of these operators: 

tr {xm (p) = ^2p v ,v 
y 

The sum converges, and the definition is independent of the choice of basis since p is 

trace class. 

Now we have the terminology to state the requirements for an axiomatic QCA. 

Definition 3.3. The global evolution R of a QCA on the Hilbert space of finite configurations 
He, with neighborhood Af, is required to be: 

(i) Unitary: R)R = I. 

(ii) Translation invariant: Let t z , for some 2; e Z", be the translation operator defined 
by its action on an element |c) = (££) JgZ n \cj) G C: 

Tz ■ (g) \Cj) H> (g) \c j+z ) 

The map t z is extended by linearity to span(C), on which it is norm preserving. 
Then t z can be unitarily extended to He ■ We denote the group of translations: 
T = {T z } z& z n - A linear operator M on He is translation-invariant if t z Mt~ x = M 
for all zeZ n . 

(iii) Causal: A linear operator M : He — > He is said to be causal relative to a neigh- 
borhood N if for every p, p' density operators, and x G Z n , that satisfy: 

p\n x = p'\m„ 

the operators MpM\Mp'M^ satisfy: 

MpM ] \ x = MpM^\ x 

We note that, just as changing the specific choice of symbols in the alphabet and local rule 
of a classical CA does not change the dynamics, any local isomorphism of the cell Hilbert 
space and the global time evolution operator of a QCA results in equivalent dynamics. The 
local isomorphisms of the cell Hilbert space are simply the unitary group of the appropriate 
dimension. It is therefore only dimension d of the cell Hilbert space that defines the local 
Hilbert space W because all Hilbert spaces of the same finite dimension d are isomorphic 
to C d . For any QCA we may therefore define a local equivalence class of QCA defined by 
a local Hilbert space of finite dimension W, all local unitary operators on W and a global 
evolution R on %c- 

Now we define a local transformation of the Hilbert space of finite configurations. Given 
a local unitary operator U acting on W whose action is: 

U : B Q — > B' Q = {U\q) : qe Q} (7) 

The set of locally transformed finite configurations Cu is defined as in Definition 12 . 1 1 with Bq 
replaced by Bq. Similarly, a local transformation of the Hilbert space of finite configurations, 
denoted by %c U i ls the I 2 completion of span(C;y), as in Definition 12. 2\ under the norm 
induced by the inner product (JSJ), which is preserved by the local unitary transformation. 

1 That such an extension exists, and is a unitary operator on Tic follows from the Bounded Linear 
Transformation (B.L.T.) Theorem, standard in the theory of operators on Banach spaces. The reader can 
find it in Reed and Simon [43] as Theorem 1.7, pg. 9. 
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Given a particular global evolution R on He an element of its equivalence class is defined 
by a transformed operator 

Ru=(0U)R(0U^ (8) 

and a locally transformed Hilbert space as defined after equation (j7j). The equivalence class 
is defined by {{R v ,U Cu ) ,VU} . 

Let us see what the constraints of Definition 13.31 tell us about an operator. First we 
consider translation- invar iance. We begin by identifying the invariants (the fixed points), in 
He, of the group of translations T = {r z } ze z n - 

Lemma 3.4. The space of T -invariants in He is one- dimensional. It is: C x &) xeZ n \qo) > 
where C x =C\{0}. 

Proof. Consider the action of T on He- Let us write <fi G %e i n a basis {|Cj) : i G Z} of He '■ 
where |C 4 ) G C. Let |C ) = ®.,. eZ n \lo)- If r z £ T, then: 

i ieZ 

Since r 2 permutes elements of C, we can associate a permutation s z of the indices « G Z with 
the action of r 2 , i.e., r z (|Cj)) = \C Sz{i) ). 

If is an invariant of T, then: r z (0) = for all z G Z n , i.e., 

£>|Ci> = ^Md)) = ^K-HolCt) 

This implies 0» = <f> Sz -i(i) for z G Z n , i.e., fa are constant on the orbits of C under the 
action of T. But T fixes \Cq) = (SWz™ \qo), and all other orbits in C are countably infinite 
in cardinality. Since <fi has a finite £ 2 norm, this implies that fa = for all i ^ 0, i.e., the 
space of T-invariants is one-dimensional: C (Q xeZn | (fo) ■ D 

Now that we know that the only vectors fixed by the group of translations are scalar 
multiples of the quiescent state: <S) xe %n \qo), we go on to determine the action of invertible 
translation-invariant operators on the quiescent state. 

Lemma 3.5. An invertible and translation-invariant operator M on He has (S^ez™ \Qo) as 
an eigenvector: 

M( (g) |g » = A (g) |g > 

xez n i6P 

for some A G C x . In particular, if M is unitary and translation-invariant, then A = e t@0 
for some Go G K. 

Proof. M is invertible, soM( <SWz" l?o)) 7^ 0- And M is translation invariant, hence: 

^^((g) |g )) =M((g) |g » V^GZ" 
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But t z 1 ((g) xeZn | go)) = <8Wz« \lo) for all z G Z n . This implies: 

t z (m( (g) |gb))J = M((g)|g » V^eZ" 

that is, M((^) ;! , gZn |go)) is an invariant of r 2 for all z G Z n . Lemma [3.41 now implies that: 

M( (g) |g » = A (g) |g ) 

for some A G C x . 

As a special case, if M is unitary then all its eigenvalues are roots of unity and so Ao = e* e ° 
for some O 6K. □ 

Finally, we consider the implication of causality on a unitary operator. First we need 
the related notion of a local operator. Consider the embedding 2 of (S^gd End(W / ) into a 
subalgebra of B(Hc) (the algebra of bounded linear operators on 1-ici where the norm is the 
usual operator norm): 



l d : (g)End(W) ^ B(H C ) (9) 
jeD 

where a is an element of <S>j &D End(W / ), and % is the identity operator on the co-D space, 
Hep, with the operator decomposition as in Through the embedding to & the algebra 
(£)j eD End(VT) is isomorphic to the corresponding finite dimensional subalgebra of B{% c ). 

Definition 3.6. An operator M on T-Lq is local upon a finite subset D G Z n if it is in the 
image of the map ld ©. 

To understand the significance of these local subalgebras, we connect the algebraic struc- 
ture of the subalgebra they generate to the topological structure of B(Hc)- Let us denote 
the subalgebra embeddings Q by A D = i D (® . D End(W)). Let {D k C Z n : \D k \ < oo} ke 



be a strictly ascending chain: 

D QD 1 C... 

such that: Z ra = U k( zf>}D k . We have an ascending chain of subalgebras formed by embeddings 
{AD k }kem, and we denote the union of these subalgebras by Z: 

Z = U k A Dk (10) 

Theorem 3.7 (Density Theorem). Z is dense in B(7-L c ) in the weak and strong operator 
topologies. 

Proof. Z acts irreducibly on "H c : it acts transitively on a dense subspace, span(C), of %c- 
This implies, by Schur's Lemma IA.lt that Comm(Z) = CI, where I is the identity oper- 
ator on He- Then Comm(Comm(Z)) = BCHc)- By the von Neumann Density theorem 
(Theorem IA.8j) . Comm(Comm(2)) is the weak and strong closure of Z. □ 



2 For any finite-dimensional vector spaces V and W, we will assume the standard isomorphisms: 
End(y) = V ® V* and End(V ® W) = End(V) ® End(VF), and make use of them as needed, without 
explicit mention. 
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The embedding i D (jUJ) allows us to consider the individual elements of A x (the special case 
when D = {x}, i.e., a single cell) as acting only on the tensor factors of interest. Then we 
are justified, for notational reasons, in replacing the finite dimensional algebras Ad with a 
finite tensor product of cell algebras of the form A x . 

For a unitary operator M on B (Jic), denote by Cm, the conjugation by M map on B(fHc)'- 

C M ■■ B{H C ) — > B{U C ) 
Z i-). M ] ZM 

As M is assumed to be unitary, this map is a unitary isomorphism (under the usual operator 
norm) of B(fHe)- 

For any set Y C B{fhi c ) the image of Y under conjugation by M is Cm(Y). By abuse 
of notation, we say, Cm(Y) = M^YM. The images of the cell algebras under conjugation 
by M, M. = .\/U .W. z G Z n are finite dimensional, and for z,z' G Z, n , z ^ z' , A4 Z , AA z i 
pairwise commute as A z , A z > pairwise commute. 

Having defined the local operators we now consider the supports of some such operators 
of particular interest. Let the reflected neighborhood of Af, denoted V, be: 

V = —Af = {—k : k G Af} 

Then this reflected neighborhood V can also be translated. We denote the translation by 
V x . 

V x = { x - k : k G AT} 

The size of this set, by definition, is |V X | = \Af\. 

It is straightforward to see, by symmetry, that V x C Z™ is the subset of the lattice 
consisting of those elements whose neighborhood contains x: 

V x = {l:xE M} 

Causality can be expressed in the Heisenberg picture in which one considers the evolution 
of operators. This form of causality is captured by the Structural Reversibility theorem due 
to Arrighi et al [3]. In the interests of making the present paper self contained we have 
included the proof of this Theorem in the Appendix (Theorem IC.ip . 

Theorem 3.8 (Structural Reversibility). Let M : Tic — > He be a unitary operator and Af 
a neighborhood. Let V = {— k\k G Af}. Then the following are equivalent: 

(i) M is causal relative to the neighborhood Af . 

(ii) For every operator A x local upon cell x, M^A X M is local upon Af x . 
(Hi) is causal relative to the neighborhood V. 

(iv) For every operator A x local upon cell x, MA X M' is local upon V x . 

The requirement of causal evolution will be useful in the local analysis of the global 
evolution R. This is because, by Theorem 13 .71 the subalgebra generated by the local algebra 
is dense in B(fhLc)- The unitary, causal evolution of the QCA can therefore be considered in 
terms of finite dimensional local subalgebras of B(fHc)- 

For a site x the set of sites whose neighborhood contains x is V x . We first show that the 
algebra A x is a subset of the algebra formed from the span of the images of the local algebras 
Ak for k G V x . This result is true for all QCAs. 
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Theorem 3.9. Let M be a unitary and causal map on He relative to some neighborhood M . 
Then for every i 6 Z", 

A x C sp&n(Il kGA fM x - k ) 

In particular Z = M^ZM. 



Proof. By Structural Reversibility, Theorem 13.81 (iii) , for every ifZ", 

MA X M^ c (g)A- fc 

This implies: 

A x = M ] MA X M^M C M f ((g) A x - k )M = span(U k€M M^A x ^ k M) = spm{U keJ ^M x . 

keAf 



The last statement follows by the first statement and Theorem 



ii, 



i.e., 



M Z C(^} A z+k for all z G Z n 

□ 

For the particular case of global evolution R of a QCA, we note the map given by conju- 
gation by R is: 

C R : B{U C ) — ► B{Uc) (11) 

Denote 7?. 2 = R^A Z R, the image of cell algebras A z under conjugation by R (fTT]) . This is 
the algebra of operators which are the images of operators localized on a single cell, after a 
single timestep of evolution by R. We are interested in the intersection of these time evolved 
algebras with the algebras localized on a single cell. Let us denote by T> Z)X the following 
subalgebra of 1Z Z , z G Z n : 

v z>x = n z nA x (12) 

When z G V x then D 2jX are the elements of 1Z Z which are contained in A x , where 

K z = RA Z R ] C End(iy) 

and: 

Ac = End(W)<8) (g) I fc . 

The subalgebras P^ j!B for z EV X may be understood intuitively as the subalgebras localized 
on x that the cell x receives from cells in V x . It turns out that the algebra generated by 
V z x for z G V x is the key ingredient in understanding the structure of a subset of QCA. 
In particular we are interested in the case when the algebra generated by T> z x for z G V x 
is equal to the full cell algebra A x . We state a preliminary theorem about the structure of 
W = C[Q] and the subalgebras T> ZtX when this is the case. This is the theorem that allows 
a local classification of a QCA as a QLGA. Let T> x ^ y ^ x be defined as in (1121) . where x G Z n , 
and y G Af. 
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Theorem 3.10. Suppose that R is the global evolution of a QCA with neighborhood M ■ 
Then A x = span(Yl y& j^V x _y jX ) if and only if there exists an isomorphism of vector spaces: 



for some vector spaces {V z } ze j\f. Under the isomorphism S , for each y G N ' : 



T^ x -y,x — End(V^) <g> Q£) I Vz 

Before proving the theorem, let us state a corollary describing the structure of the algebras 
7l x under the conditions of the theorem. These algebras TZ X are the images of algebras A x 
localized on a single cell, after one timestep of the global evolution R. 

Corollary 3.11. Suppose that R is the global evolution of a QCA with neighborhood Af , and 
satisfies A x = span(Jl ye j^V x _ yx ) ■ Then: 

(i) A x = End(W) (g) 2gA ,End(X) ; for all x G Z n . 

(ii) The dimension of W, dw, is a product of the dimensions of V z , dy z , i.e., d\y = 

(iii) 7l x = span(U k£Af V X:X+k ) £ (g) A;gA/ -End(V fc ). 
Proof. 



of 



in 



and (ii) are obvious from the isomorphism S in Theorem 13.101 For the proof 
we note that by definition of T> x ^ x+k (IT2]) . 1Z X D spa.n(Jl k€ j^V x , x +k)- But from 
Theorem EDJ span(U keM V x ^ x+k ) ^ <g) k( ^End(V k ) = End(V^). We already know that 
TZ X = R)A X R = End(W). Hence the conclusion. □ 

This corollary shows that if the conditions of Theorem 13.101 are met then the algebra A x 
is isomorphic to an algebra generated by a set of commuting algebras each of which is the 
endomorphisms (full matrix algebra) of a single tensor factor of the space W. 

The following lemma shows that if there is a set of pairwise commuting self-adjoint algebras 
that generate End(W), then there is a tensor product decomposition of the space W such 
that each of the commuting self-adjoint algebras is the set of endomorphisms (full matrix 
algebra) of one of the tensor factors. This will be useful in the proof of Theorem 13.101 

Lemma 3.12. Let {M k : M k C End(W / )}^ =1; be a finite set of distinct pairwise commut- 
ing self-adjoint algebras, each containing the identity operator, such that span(n^ 1 Mfc) = 
End(W / ). Then there is a vector space isomorphism: 

m 
k=l 

for some vector spaces {Vfc}^L 1 . Under the isomorphism Sw, for each 1 < y < m: 



M y = End(V y ) (8) (X) l Vk 

fce{i,...,m}\M 

Proof. If m = 1, there is nothing to prove. So assume m > 1. Consider M m from the set 
{M k : M k C End(H/)}^ =1 . Let T m = span(n[^ 1 Mt). Then we know that T m and M m 
are mutually commuting self-adjoint subalgebras, each containing the identity operator. By 
Proposition IB.2} W is then a completely reducible T m -module. Corollary IA.6I then implies 
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that there exists a finite set of irreducible T m -modules 3 that we denote {VjK =1 , such that 
if we denote the set of their multiplicity spaces: {Uj = Hom Tm (V ? -, W)} r j =1 , then we have a 
T m -module isomorphism S m : 

r 

S m :0f/ 3 ® Vj — »• W 

3=1 

r r 

J2 U 3 ® V 3 ^ ^2 U A V l) 

3=1 3=1 

Under this isomorphism: 

r 

T m = S m ( l v . <g> End(V^)) S' 1 

3=1 

r 

Comm(T m ) = S m { End(C^) <g> ly^S' 1 

3=1 

where Comm(T m ) is the commutant of T m , as defined by (HT1) . and each Uj is a Comm(T m )- 
module under the action given by fj42l . Since, by definition of T m , M m C Comm(T m ), this 
implies: 

r 

span(M m T m ) C span(Comm(T m )T m ) = S m ( End(V^) <8) End(t/ i ))5 m 1 C End(W) 

i=i 

But by assumption, span(M m T m ) = End(W / ). Then the inclusions above are equalities, and 
there must only be one summand in the sum above. That is, there exists some irreducible 
T m -module Y m , such that if we denote its multiplicity space: U m = Homr m (F m , W), then we 
have a T m -module isomorphism S m : 

S m :U m ®Y m -^W (13) 

Furthermore, under 5,^: 

T m ^ I Um (g) End(F m ) (14) 

and: 

M m = End(C/ w ) (g) Iy m 

Here the subscript m signifies that Y m is the only class of irreducible T m -module in W, 
and does not correspond to the indexing {j : 1 < j < r} used above. This subscripting 
convention will be used for general k in the set: {k : m — 1 > k > 2}, in the arguments to 
follow. 

In particular this proves the lemma for m = 2: if m = 2, observing that T2 = Mi, we can 
define the vector spaces, V2 := U2, and V\ := I2, and the isomorphism Sw '■= S 2 X - Then the 
vector spaces V± and V2, and the isomorphism Sw satisfy the statement of the lemma. 

Now assume m > 2. Since T m is a subalgebra of End(W / ), we choose the irreducible 
T m -module Y m in (fT3|) to be a subspace of PU: F m C W. 

Clearly the subalgebras in the set {Mfc}^ 1 are subalgebras of T m = span(n^ 1 M t ). 
Also, we have that T m = End(F m ), by the definition of the T m action on Y m , and (JH]). 

3 Here the index set { j : 1 < j < r} enumerates the finite (by Corollary IA.6|) set of equivalence classes of 
irreducible representations of T mi denoted T m . Each value of j corresponds to a distinct class A G T m . 
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Let us consider the restricted actions of the subalgebras in the set {M^}^ 1 on Y m . These 
restricted actions comprise a set of mutually commuting, self-adjoint (for the inner product 
on W restricted to Y m ) subalgebras of T m = End(y m ), each of which contains the identity 
operator on Y m . 

We proceed by induction for this part of the argument. Let k be in the set: m — 1 > k > 3. 
LetT fc = span(n*LiM t ). By way of induction, assume that there exists a descending sequence 
of subspaces of Y m , denoted {Y k }^Z^ satisfying: Y m D Y m _ x D . . . D Y 3 , such that each Y k 
is an irreducible T^-module under the restricted T k action. Since T m D . . . D T 3 is a 
descending sequence of subalgebras, therefore, T k can be restricted to act on Y k+ i. Let 
U k = B.om.T k {Yk,Yk + i). Also assume that there is a T^-module isomorphism: 

Sk :U k ®Y k ^Y k+1 (15) 
u k ®y k ^>- u k {y k ) 

where u k G U k , and y k &Y k . Under this isomorphism, assume that: 

T k = I Uk ® End(n) (16) 

and: 

M k End{U k ) ® l Yk 

Now let us consider the terminating case when k = 2. The subalgebras Mi, M 2 are sub- 
algebras of T 3 = span(MiM 2 ). By induction assumption ffl6l) . under restricted T 3 action, 
T 3 = End(Y"s). Restricted actions of the subalgebras Mi, M2 comprise a set of mutu- 
ally commuting, self-adjoint (for the inner product on W restricted to I3) subalgebras of 
T3 = End(Y3), and each of them contains the identity operator on Y3. Also, T2 = Mi. 
This case then reduces to the case m = 2, that we have already proved before stating the 
induction assumption. The m = 2 result then implies that for the restricted Mi action on 
Y 3 , there is some irreducible Mx-module Y 2 , such that if we denote its multiplicity space: 
U 2 = B.om Ml {Y 2 , Y3), we have an Mx-module isomorphism S 2 : 

S 2 :U 2 ®Y 2 ^ Y 3 (17) 
u 2 ®y 2 ^ u 2 (y 2 ) 

where u 2 6 U 2 , and y 2 e Y 2 . Under the isomorphism S 2 , we have that: 

Mi = I U2 ® End(F 2 ) 

and: 

M 2 = End(^ 2 ) ® Iy 2 

Using the isomorphisms (fT5l) of the induction step, with (TlBl . and (fTTj) . implies that S m 
is given as: 

S m :U m ® E7 m _! ® • • ■ <g) ?7 2 g) r 2 — >• 

u m <S> u m -i <g> ■ ■ • <g> m 2 ® yi ^ w m (w m _i(. . . (u 2 (y 2 )))) 

where Mfc G ^4, for 2 < A; < m, and y 2 G y 2 . 
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Define the vector spaces V k := U^, for 2 < k < m, and V\ : = Y 2 , and the isomorphism 
Sw '■= S' 1 . Then the vector spaces {Vk}™ =1 , and the isomorphism Sw satisfy the statement 
of the lemma. 4 ' 5 

□ 

Remark. We observe that if some of the subalgebras M k in Lemma 13.121 are multiples of the 
identity operator, then the corresponding vector spaces V k will be (trivial) one-dimensional 
spaces. 

Proof of Theorem \3.1CK Let us assume that A x = span(Jl ye j^V x _ yx ). Fix x G Z n , let y G Af. 
The subalgebras T> x _ yx are self-adjoint (as lZ x - y are self-adjoint) subalgebras of A x , and they 
pairwise commute (as TZ x - y pairwise commute), and each contains the identity of A x . Hence, 
by Lemma 13.12} there exists a isomorphism of vector spaces: 

S x : W-+®V x - k , x 

keAf 

for some vector spaces {V x -k, x }k&N- Under the isomorphism S x , for all k G Af: 

V x - k , x ^End(V x „ k>x )® (g) l Vx _ z>x 

z£Af,z^k 

By translation invariance of R, the isomorphisms S x can be taken to be the same and the 
sets of vector spaces {V x _k,x}keN, f° r an x £ can De taken to be the same. Therefore, 
there exists a isomorphism of vector spaces, S: 

for some vector spaces {V z } z£ ^/. Under the isomorphism S, for all y G Af: 

V x _ y , x ^End(V y )® (g) % 

z£Af,z^y 

For the converse, let us assume that there exists, for some vector spaces {V z } z& _\f, an 
isomorphism S: 

S:W^0V z 

zeAf 

such that 

V x - y , x End(y„) ® (g) Iv, 

zeAf,z^y 

then: 

span( J] £> x _ y J = S^spanf J] (End(V„) ® (g) lyJ^S 

j/eA/" ^yeAf zeAf,z^y ' 

= S- 1 {g)End(V z ))S = A x 

zeM 

4 The proof shows that the order in which the subalgebras {Mfc}£L;L are enumerated does not affect the 
conclusion. 

5 One can also prove this lemma via a slightly more general route, noting that the subalgebras involved are 
semisimple, which follows by Corollary IB. 3[ and by using the semisimple version of the Double Commutant 
Theorem, Theorem I A. 71 Approach in the given proof is more specific. 
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□ 

Theorem 13.101 establishes an equivalence between the fact that the algebra of operators 
on a single cell is generated by those pieces of evolved subalgebras that are localized on the 
cell, and a tensor product structure of the single cell Hilbert space. In the following we show 
that Theorem 13.101 enables the classification of a subset of QCAs as QLGAs. 

3.1. Separability of the quiescent symbol and the propagation operator. Assume 
there is an isomorphism S as in Theorem 13.101 Denote W = S(W) = <S> zeJ ^V z . W can 
be identified with C dw and imbued with an inner product. We choose this inner product 
such that the image of the symbol basis Bq (TJJ under S, S(Bq) = {S(\q)) : q G Q} is an 
orthonormal set. With this inner product, W becomes a Hilbert space. Note that this basis 
of W, which is the image of Bq, will not in general consist of simple tensors. In addition, 
imbuing each Vk with an inner product and then constructing an inner product on W may 
result in an inner product on W which is inconsistent with the requirement that the image 
of the symbol basis be orthonormal. 
For each z G Af, let: 

{\k z )}i< kz < dVz (18) 
be a basis of V z . Use k to label elements of the Cartesian product U ze _\f{l, .., dy z }- A natural 
basis for the Hilbert space W = <S> ze j^ V z can be written: 

B N = {(g) |k(*)) : k G IW{1, dv z }} (19) 

The notation k(z) stands for the 2-th coordinate element of k G U z£ n{1, .., dy z }- And where: 

|k(z)) G V z (20) 

We shall refer to an individual tensor factor within a cell, i.e., an element of V z for some 
z G A/", as a component. The notation here is just as for the case of a classical LGA in 
equation (JTJ). In the classical case different lattice vectors z can have different alphabets. 
Here, the dimensions of the V z may differ, and hence they have distinct bases. In the case 
where V z = C 2 V z the basis of W that we obtain corresponds to the occupation number 
representation for systems of particles. This representation is used widely, particularly in 
quantum simulation of fermionic systems [5,47,59]. 
We add \q ) = S(\q )) to the basis B^ and define: 

Q = B Af U{\q )}. (21) 

This is necessary because \qo) may not be an element of the basis B, in which case we can 
express it as: 

\q )= "(k) <g) |k(z)> 

k 1 1 .\ ! • •••. , '\ } zeAf 

where k G n z& y{l, .., dv z } and u is some complex valued function: 

u : n 26Ar {l, ..,d v J — ► C 

We then consider the set of finite configurations, taken over every site, and define the 
counterparts of the set of finite configurations and the Hilbert space of finite configurations: 

Definition 3.13. The set of finite configurations in component form, denoted by C, is the 
set of simple tensor products {(£) xeIi n \c x ) : \c x ) G Q, all but finite \c x ) = ]%)}. 
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The inner product on C is induced by the inner product on W. Let |c) = (S>xez™ l£c)> 1^') = 
<S> xeZ n \c' x ) G C. Define the inner product of the elements \c), \c')\ 

(c|c') = U xeZ n(c x \c' x ) 

and extend it by linearity to get an inner product on span(C). 

Definition 3.14. The Hilbert space of finite configurations in component form, denoted by 
% c -, is the £ 2 completion of span(C), under the norm induced by the above inner product. 

We will use the term Hilbert space of an individual cell x to simply mean the set in 
which | Qr) can take any value, without violating the definition of C. The tensor factors 
corresponding to the Hilbert spaces of individual cells are indexed by the cell index x G Z n . 
Thus we can index the Hilbert space of a cell x by W x := W . These spaces themselves 
are composed of tensor factors V z indexed by the component index z, which runs over the 
neighborhood Af, i.e., z G Af. So we index the sub-factors of W x by V ZjX := V z . Then 

w x = ® zeAf v g>x . 

The picture of time evolution that we shall establish is that of transformations acting 
only on each cell Hilbert space ^ ze j^V z , followed by a permutation on tensor factors that 
rearranges the components Ik^z)). We first take care to define the map that permutes the 
tensor factors and understand what its existence implies about the structure of the space 
%£. These properties will be useful in proving the main theorem. 

We will formally specify the propagation operator, relative to the neighborhood A/", by 
specifying it on elements of C. This map acts by shifting one component from each of the 
the basis elements (£) ze j^\k y (z)) of the neighborhood cells y G Af x into the corresponding 
component of the cell x. 

We see at once that to be able to define the propagation operator on C will require that 
the new quiescent symbol \q ) = S(\q )) G (£) ze xV z , is a simple tensor, i.e., of the form: 

l?o) = !*>.*) 

z£jV 

for some \%,z) G V z . If this is not the case then we may always find a unitary operator U 
acting on W, such that S(U\qo)) is a simple tensor. In this case we work with the locally 
equivalent QCA defined by the pair (Ru,7icu) in (JED- Let us therefore assume that |g ) has 
this property. 

Under this assumption, the basis for V z in (Ti~8]) can be chosen such that \qo) G B^f. With 
the inclusion of \q ) in Bjy, the definition of Q (12"T|) becomes Q = (IT91) . This allows us 
to write an element |c) G C in the form: 

|c) = %) =00 M*)) ( 22 ) 

xez n x<=z n zeAf 

where the arth tensor factor of an element (£) xeZ n \c x ) G C is given by: 

\c x ) = |k,(^)) G W x (23) 

z&N 



and where: 



K(z)) G v Zy 



(24) 
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Now we may define the propagation operator relative to the neighborhood Af, that we 
denote by a. We combine the cell and component indices together in a pair (x, z) G Z n x N '. 
As a permutes the cell indices x G Z n , leaving the component indices z G M intact, we can 
define a, on C, in terms of a cell-indexing function s a : Z n x M — y Z 71 , as: 

a : C > C / 2 ^n 

where s a is defined as follows. 

s a : Z n x AT — ► Z n 

(x, 2) X + Z 

In other words a moves the z-th component of cell x + z to the z-th component of cell x. 
a is extended to a unitary, translation-invariant operator on H c . We can omit s a and write 
the action of a as: 

o- : ®* e z» <8W l k *0)> ^ <8>* e z« <S), e ^ 
We note that the structure described above corresponds to the classical case in the follow- 
ing way. We identify the states of a classical lattice gas with the elements of C. The state 
of a lattice gas site is simply \c x ), and the state of lattice vector z is an element of the basis 
of V z . We may think, therefore, of the tensor factors V z , as Hilbert spaces of lattice vec- 
tors pointing to neighborhood site z. Note however that because these are quantum models 
arbitrary superposition states over the elements of this basis are allowed. 

3.2. The collision operator. As we shall prove below, an evolution composed of the cell- 
wise isomorphism S and the propagation operator cr obeys the conditions to be a QCA. 
However, this is not the only possible evolution. Performing any cell-wise unitary after the 
isomorphism S will also give us a QCA. 

We therefore formally construct an operator F, by specifying it on the basis C of in 
terms of some unitary operator F on W = <S> z( zj^V z , as follows: 



F : ® |cA ^ (g) F{%)) 



Then we have the following lemma that tells us what conditions ensure that F is defined on 
Ug: 

Lemma 3.15. F is defined on C as a map F : C — > H^, and can be extended to as 
unitary and translation-invariant operator, if and only if F has \%) as an invariant (an 
eigenvector with eigenvalue one): F\q ) = \q ). 

The proof of this lemma is included in Appendix [Dl Suppose F\q ) = \q ), then we can 
define: 

F : — > He 

which is the extension of the following definition on C to Hg. 

F:C—^H C (26) 
%} ^ 6d F(\c x )) 



x£Z n 
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Whereas the operator a is a map that moves the tensor factors V z between cells of the 
automata, the operator F acts locally on states in <S> zejK fV z . In the previous subsection we 
identified the tensor factors V z with lattice-gas vectors. The operator F acts separately on 
the states of each lattice gas site. This operator corresponds to the collision operator of a 
quantum lattice gas, as we shall show in detail below. 

3.3. The definition of a Quantum Lattice Gas Automaton (QLGA). At this point 
we have developed tools that are sufficient to give a formal definition of a Quantum Lattice 
Gas (QLGA). The dynamics proceed by collision and propagation, as in the classical case 
and as in the models defined in [34,36-38] and [8-10]. A QLGA is defined as follows: 

(i) A lattice C = Z n , whose sites we shall label by x. 

(ii) A neighborhood M . This translates to a neighborhood M x for each cell x, consisting 
of a set of cells as follows: Af x = {x + z\z G Af}. 

(iii) For each element of the lattice £ there is a site Hilbert space W = ^ zeJ ^ V z , for some 
finite-dimensional vector spaces {V z } z ^. W contains a distinguished quiescent state 
| go) which is a simple tensor: 



(iv) A Hilbert space of finite configurations in component form Hq as defined in Defini- 
tion EH 

(v) A state of the automaton, is an element of the Hilbert space of finite configurations 
in component form. 

(vi) A propagation operator relative to the neighborhood A" as in ff25l) . 

(vii) A local collision operator F, which is a unitary operator on the site Hilbert space 
W = <S> ze j^V z , such that F has |go) as an invariant (an eigenvector with eigenvalue 
one): F\%) = \%). 

(viii) A global evolution operator R which consists of applying the propagation operator 
o followed by the local collision operator at every site, i.e., is given by: 



where F is the collision operator constructed from F as in f )26|) . 

The tensor factors V z may be identified with the lattice vectors of a classical lattice gas, 
which point to the elements of the neighborhood of a site. 

3.4. Local inner product preserving transformations to component form. Let S 

be the map that takes He to "H^, i.e. the global map that is the result of applying the local 
isomorphism S to every cell. 



Qo) = (g) \%,z), where \q 0>z ) G V z . 



Z&JV 



R = Fa 



S:H C 




this map is defined on C as: 




and extended to He- Here S : W — > W is as in Theorem 13.101 and with the choice of inner 
product onW = <S> zeJ ^ V z as discussed at the beginning of Section I3"TT| S is an inner product 
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preserving isomorphism of Hilbert spaces W and W. This implies S is an inner product 
preserving isomorphism of Hilbert spaces Tic and Tig. 

3.5. Quantum cellular automata and quantum lattice gases. Now we combine the 
foregoing results to prove our main theorem and characterize the class of QCAs that are 
QLGAs. For any QCA defined by (R, Tig) we can construct a locally equivalent QCA 
(R, Tic) that satisfies 



l<?o) = S(\q )) G QQV Z , is a simple tensor, i.e., \q ) = Q£) \q 0jZ ), where \q 0)Z ) G V z . 
zeM z^M 

We therefore consider a QCA (R,Tic), with neighborhood Af, and for which |g ) is a simple 
tensor. Let T> x _ yx be defined as in (1T2)) . where igZ", and y G Af. 

Theorem 3.16. R is the global evolution of a QCA on the Hilbert space of finite configura- 
tions Tic, with neighborhood M , and satisfies: A x = span(U ye j^'D x -.y )X ) , if and only if: 

(i) There exists an isomorphism of vector spaces S: 

S:W — >(g)V z 

zeN 

for some vector spaces {V z } z& j^. Under the isomorphism S , for each y G M : 

V x ^ x = ^{Vy)® l Vz 

Furthermore, W = ® zeJ \f V z can be given an inner product such that S is an inner 
product preserving isomorphism of Hilbert spaces. 

(ii) R is given by: 

R = §- l Fa§ 

where a, F, and S are as in (122]), ( 125]) ■ and ( 127)) respectively. 
(Hi) F, in the definition of F, is unitary and has \%) as an invariant: F\q ) = \q ). 

Proof. We first prove that if R is the global evolution of a QCA, and satisfies: 

A x = span(U ye j^V x -y tX ), 



then (i) , (ii) and (iii) follow. 



By Theorem 13.10) the vector space isomorphism S in (i) is immediate. We can imbue 
W = S(W) = (£) ze j^V z with an inner product such that {S(\q)) : q G Q} is an orthonormal 
set (as discussed in the beginning of Section I3.ip . By this choice of inner product, W is 
a Hilbert space such that S is an inner product preserving isomorphism of Hilbert spaces, 
hence proving ( 



By corollary 13.11 ((iii) TZ X = span(U z( zj^V x tX + z ), where the definition of T> x>x+Z is as in (j!2p . 
It follows that: 

STZxS' 1 = span(n ; , gA ^S'P XiX+;Z S' _1 ) (28) 

where S is as in ( 12 7p . We use the indexing from ( I22p . (123 p . and 023]). The isomorphisms 
from Theorem 13.101 give us: 

A x = S- 1 (0End(V z , x ))S 

z&N 



SV^S- 1 = End(V 2 



Z,X + Z ) } 
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Using the above, we expand the LHS of (1251) : 

STZj- 1 = SrfAxRS- 1 = {Srf End^)) {SRS' 1 ) 

z<=Af 

and the RHS: 

span(n^ eA r£'D :E>:c+2 S , ~ 1 ) = span(U zeA fEnd(V z , x+z )) = (g) End(V Z)X+z ) 

and obtain: 

(SrfS- 1 ) ( (g) End(K,s)) (.mS- 1 ) = (g) End(K,*+*) (29) 
zeftf zeAf 

Let us denote R = SRS^ 1 . 

R — S RS : Hq — y ^~^C' 

R is clearly unitary as R is unitary and S is inner product preserving (S is composed of 
inner product preserving operator S acting on every cell), and is translation- invariant since 
R and S are translation-invariant. 

Fix x G Z n . Choose a basis, By x , of W x = &) zeA f V z>x : 

4 = Mcg^ (30) 

z£jV 

For instance, we can choose B Vx to be the product basis of W x of the form (JHJ), letting : 
= {® zeA f l k *0)> ■ k - e n zeAr {l, .., d Vz }}. The set ® (*/J : G B % } is a 

basis of End(M4) = (§) zeA/ - End(V^ )iZ ;). Choose such a basis element of End(Wx) and create a 
basis operator local on cell x: 

A x = {\b x )®{b' x \)®l-^ (31) 

where Ij^t is the identity on co-{a;} space (defined in Definition 13.11 for He, and the definition 

applies equally for Tig by replacing C with C). It^t is included as part of the local operator 
to help make the argument clear . 

R}A X R = R)({\b x ) ® (b' x \) ® I M )R 

By ( 1291) . i^Agi? is a rank one element of (^) 2gA ^End(\4 iX+2 ), whereas, by definition, A x 
is a rank-one element of End(W / x ) = ® ze ^f Fmd(V Z)X ). Since WA X R is an element of 
(3) zeAr End(V ZiX+2 ) for all basis elements A x G End(W / x ) of the form chosen, this implies 
that we can write W A X R as follows: 

R)A X R = 4({\w x ) ® «|) ® %)<r x . (32) 

for some G = (^) zgjV -K,x, and an appropriate permutation operator a x . o x 

maps the components of the neighborhood elements of cell x, <£) zeJ ^V x+Z)X , to the corre- 
sponding components of cell x, <S> zeJ ^ V Z;X . 

We label the combined cell and component indices together by the pair (t, z) G Z n x N '. 
It is necessary and sufficient that a x satisfy the following conditions: 

(a) a x permutes the z-th component of a cell to 2-th component of another cell (where 
z G AT). That is, a x permutes the cell indices t G Z n , leaving the component indices 
z G Af unchanged. 



WHEN IS A QCA A QLGA? 



25 



(b) For the cell index t = x, a x maps the z-th component of the neighborhood cell x + z 
to the z-th. component of cell x (where z G AT). 
To satisfy the requirement in condition (a) that a x only permute the cell indices, keep- 
ing the component indices intact, we use a cell-indexing function, s x : U 1 x M — > Z n , 
in the definition of a x . As every element \c) G C (Definition 13. 13j) is of the form |c) = 
<S>tez n ®zeW l^t^)) ( l22|) - we can define a x as follows: 

a x : C — > C 

®t G z"8WM*)> ^ <8) ieZ n0 ze Ar \K(t,z)( z )) 

The rest of the requirements in conditions 



s x . The remaining requirement in condition 



and (b) above can be recast as restrictions on 
'a) is that o x is a permutation. This implies 
that when the component index of the domain of s x is fixed to be some z G A/", s x is a 
bijection of the cell indices Z n , i.e, a permutation of cell indices. 



s x 



Z n x z< — ► Z n 



Condition (b) implies that when the cell index t G Z n of the domain of s x is fixed to be x, 



s x is as follows: 



Sx \xxAf- X X Af > Z™ 

(x, z) i-> x + z 

a x thus defined on C under the constraints 



and (b) , is extended to "H^. 



Now consider the implication of combining ( 13T|) and ( 132|) : 

#((!&*) ® (O = 4((k*> ® («4I) ® %K (33) 

We have the above equality for every \b x ) g) (6 X | and for the corresponding \w x ) <g) (u^l- So 
we define a linear map on = <S> z eAf ^,x, by defining it on the basis Bv x ( J30l) : 

: — ► 

\b x ) i ^ |w x ) 

where and jm^) are as in (133]) . is extended to by linearity. Then T x is unitary 
since i? is unitary. By symmetry, we can write ( 133|) as: 

£ f ((K> ® (O ® = ® W2) ® %H 

This is true for every x G Z n . In addition, i? = SRS~ 1 is translation- invariant, which implies 
that all the maps T^. can be taken to be the same unitary map T on W = <S) z eM anc ^ that 
the map cr^ for each x is translation- invariant, i.e., each a x is the map a defined in fl25l) . Now 
define F : \— > W, to be F := TL Then i? = SRS^ 1 is a composition of the propagation 
operator a followed by the collision operator F: 

SRS^ 1 = Fa 

where a is as defined in (|25|) . and F is a unitary map on as defined in ([26]) in terms of 
F. Thus we get part (ii) of the theorem. By Lemma 13.151 we get part (iii) of the theorem: 
F in the definition of F has \cjo) as an invariant, i.e., F\qo) = \qo). 



Next we prove that if 



(i). 


(ii) 


and 


(iii) 



neighborhood A/", and that it satisfies the condition A x = span(n yg _/^P x _ ?/j:r ). Unitarity and 
translation-invariance are obvious from the definition. To prove causality, let us assume A x 



2(> 
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is a local operator on cell x. We consider the action of W A X R on basis elements of He- Both 
S and F act locally on cells by S and F respectively The conjugation of A x by S is an 
isomorphism of A x to A x = <^ zeJ ^End(V ZjX ) (the counterpart of A x ), and the conjugation 
of A x by F is an isomorphism of A x to itself. Observe also that F\q ) = \q ). Then we may 
as well simply look at the conjugation of A x by a, and only need to show that propagation 
is causal. Therefore we can simply consider the action of a A x a on basis elements of T-Lg. 
We can write such a basis element as: 

I3 = <8>«8>M*)» 

y&L n zeAf 

Then: 

^> = (g)((g)IV^)» 

There is no loss of generality (by linearity) in assuming that A x G A x is a basis element: 

A x = (((g) Mr))) ®(®«(r')|)) g)I M 

\ rGA/" r'eJV ' 

where (^(r)) G and \a' x (r')) G V r /, and where Ij^r is the identity on co-{x} space, as in 
Definition 13.11 Denote by p x = H ze ^(a x (z)\k x+Z (z)) . Then this implies: 

a- 1 A x a\c)= Px (g) (\a x {z))® (g) 1^'))) <g> ( (g) (g) |k y (^')» 

y=x+z-.zeN z'eN\{z} yez n \{x+z-.zeM} z'eM 

The above shows that a~ l A x a is local upon J\f x , hence by Theorem l3.8l (Structural Reversibil- 
ity), R is causal relative to neighborhood M . This in turn implies, by cell- wise conjugation 
with local isomorphisms S and F, that R is causal relative to neighborhood M . 

Next we show that A x = span(Il 3/6 _^-'D 2 ._ 3/ -,,). By the same reasoning used in showing 

causality, we can consider the following quantities. Let B z = o~ x A z o (the counterpart of 
1Z Z ). Let V x - y , x = B x -y H A x (the counterpart of U x _ ytX ), for y G Af. Then it is clear by the 
definition of a, that: 

A- y ,x = B x -y n A x = End(V y>x ) ® (g) ly^ 

z£j\f,zjty 

for each y G N '. This implies: 

X = span(L T j /6A ^P a; _j /ia; ) 

This further implies, by cell-wise conjugation with local isomorphisms S and F, that: A x = 
span(Il yeA fV x _ y)X ). □ 

The significance of this theorem is that it connects a local condition to the global structure 
of the QCA. We state the theorem in terms of propagation and collision operators which are 
the quantum analogues of the substeps of the evolution for a classical lattice gas. There are, 
however, a few equivalent ways of expressing the global evolution R under the conditions of 
Theorem 13.161 

Corollary 3.17. The global evolution in Theorem \3.16\ can be given as: 

R = faS 



for some map T: 



Proof. Let: 
on C as: 
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T : Hq — y He 
(g) \c x ) ^ (g) s-^dO) 
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S, and F are as in Theorem 13.161 T is extended to H^- 
By Theorem 13 . 161 : 

R = faS 

□ 

This is called the two- layered brick structure [45]. However, instead of working with the 
space He, we can look at the QCA global evolution as it acts on the space %£. This gives 
us another equivalent form for the global evolution R. 

Corollary 3.18. The global evolution in Theorem \3.16i can equivalently be given by a unitary 
operator R on : 

R = Fa 

Proof. Let: 

R : 7~Lq — y 

R = SRS- 1 

By Theorem 13. 16t 

R = Fa 

□ 

Thus the global evolution R is equivalent to one given by another global evolution oper- 
ator R, with S and S^ 1 providing the change of basis between He and Hq. The equivalent 
evolution R happens in two stages: the propagation (or advection) stage given by a, followed 
by the the collision (or scattering) stage given by the map F on <S> zeA fV z . This completes 
the proof of equivalence of QCA satisfying the local condition of Theorem 13.161 to Quantum 
Lattice Gas Automata (QLGA). Theorem 13.161 therefore characterizes the QLGA as a sub- 
class of QCA. It is clear that the evolution given by Theorem 13. 161 has as special cases both 
the classical lattice gas defined in the introduction and the lattice gas models previously in- 
vestigated [8-10,34-39]. In the following section we will give some examples of QCA which 
are not lattice gases and show how they violate the conditions of Theorem 13.161 We also 
compute the matrices T> Xty explicitly for the simplest quantum lattice gas and show that the 
condition is indeed satisfied. 
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4. Examples 

In this section we consider examples of three subclasses of QCA. In the first case we 
examine the simplest quantum lattice gas, defined as an example in [17] and also obtained 
from consideration of the absence of scalar QCA in [34] and investigated in [34,36-38]. 
Here it is instructive to see the algebras that appear in our local condition. Secondly, we 
investigate the counterexample of [3] to show that our condition is indeed not satisfied in 
this case. Finally, we give an example in which quantization of a classical automata changes 
the neighborhood size, and results in a QCA which is also not a QLGA. In this case we show 
how adding a propagation step to the original QCA allows our condition to be satisfied and 
hence results in a QLGA. 

4.1. The simplest QLGA. The simplest QLGA was originally investigated by Meyer as an 
extension of a single particle partitioned QCA [34]. In this model each site has two neighbors, 
and each site has two lattice vectors that point to these neighbors. At most one particle 
may be present on each vector so the site Hilbert space W is isomorphic to C 4 ~ C 2 <g C 2 . 
The collision operator F is particle number conserving - meaning that it is block diagonal 
with two blocks of dimension 1 and one block of dimension 2. We may choose the basis 
{|00}, 1 1 } , 1 10} , |11)} for W in which the collision operator F may be written: 

/l 

_ ie ia sin 9 
e ia cos9 
\0 

The lattice is Z and the neighborhood is M — { — 1, +1}. Because we are beginning from 
a QLGA dynamics composed of propagation and collision we may take the map S to be the 
identity. Then 

W = W = C 2 <g> C 2 (35) 

and the two tensor factors are V_i ~ C 2 and V\ ~ C 2 . From the Heisenberg point of view 
of the dynamics of the algebras A x , which are localized on a lattice-gas site and hence act 
on C 4 , the collision step is simply some isomorphism of the algebra. We need therefore only 
to consider the action of the propagation operator to determine the algebras TZ X and T> XjV . 
The action of a on an operator localized on the site may be written: 

a' 1 ! <g> B <g> I <g> I <g> A <g> la = I <g> I <g> A <g> B <g> I <g> I (36) 

where the first two tensor factors (reading from left to right) are the left lattice-gas site 
x — 1, the middle two tensor factors are the middle lattice- gas site x and the right two tensor 
factors are the right lattice-gas site x + 1. From this it is straightforward to write: 

£>_ M = I® End(C 2 ) 
V +hx = End(C 2 ) (g I 

Evidently the span of the product of these algebras is End(C 4 ) and hence this model satisfies 
the condition to be a quantum lattice gas, as it must. 





cos 9 
* sin 9 




0\ 




pi/3 



(34) 
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4.2. QCA that are not QLGA. Consider a one- dimensional classical CA as in Figure H] 
in which a cell is 3 bits. Then a cell x is given by a triple (a x , b x ,c x ). The neighborhood is 
Af = { — 1, 0, 1}, and the update rule \i: 

/i : o,4 a x + b x - X 
b x b x 
c x i— y c x + fe^+i 

Quantizing this CA, we get a QCA. That means, we take the rule above to describe the 
evolution of a cell state as a function of neighborhood state. A cell x of the QCA consists of 3 
qubits: \a x ) Cg> \b x ) £g> \c x ), where a x , b x , c x G Z/ (2). The quiescent state is \q ) = |0) eg) |0) <g) |0). 
The neighborhood is still Af = { — 1,0, 1}. 



4 i. 



ai 



&i 



c-i 



a 2 



4^ 



4 



r it 



C-2 



a 3 



c 3 



Figure 4. One-dimensional CA which is a not a QLGA after quantization 



Define the global evolution R on He through its action by a controlled-NOT operation on 
the neighborhood cells as follows. 

R : \a x ) \a x + 

l&x) >-»■ 16*) 

|cx) ^ |c x + b x+1 ) 

We can equivalently look at the qubits in state: |(— l) ax ) ® \b x ) <S> \(—l) Cx ), where a x , b x , c x G 
Z/(2). Since |+) and |— ) are characters for the shift representation, R can be given as: 

R: \(-l) ax ) ^ \{~l) ax ) 

\b x ) ^ {-l) hx{ax+1+Cx - l] \b x ) 

K-i)-)^l(-i)-) 

We now compute the subalgebras U x _ y)X = 1Z x - y D A x , where y G Af. V XtX = 1Z X fl A x = 
span({|a :E )(a a; | ® \b x )(b x \ ® \c x )(c x \}), where | 

) j I C X ) e {|+),|-)}, and 16,) G {|0),|1». 
£> X ._ M = TZa-i n Ax = CI ® I ® I. D x+ i,, = 72^+1 nA = Cl8l®I. Then it is clear that 
dim span(ilj /e {_ li o i i}7^x-j/,x) = 8, while dim A x = 64. Theorem 13.161 then implies that this 
QCA is not a QLGA. 

The 2-dimensional version of the above example is referred to as the Kari CA in Arrighi 
et al [3] . This is illustrated in Figure |5j The lattice C = Z x Z. A cell of this classical CA 
consists of 9 bits. The neighborhood is Af = { — 1,0, 1} x {—1,0, 1}. The center bit of the 
center cell acts by controlled-NOT on the appropriate bit of the neighboring cell as shown. 
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Figure 5. Two-dimensional Kari CA which is a not a QLGA after quantization 

In the quantized version, the cells of the QCA are 9-qubits each. Enumerate the qubits 
of cell x by (i x ,j x ) G { — 1,0,1} x { — 1,0,1} (not to be confused with the neighborhood 
coordinates on the lattice). (i x ,jx) — (0,0) corresponds to the center qubit of a cell. The 
quiescent state is |go) = ^(! l3l )e{-ioi}x{-ioi} l^)- The center qubit of the center cell acts 
by controlled-NOT on the appropriate qubit of the neighboring cell as for the classical CA. 
This action is exactly as in the 1-dimensional case. The neighborhood is the same as for the 
classical CA, Af = {-1, 0, 1} x {-1, 0, 1}. 

We compute, as in the 1-dimensional case, the subalgebras T> x _ yx = lZ x - y fl A x , where 

y G N ■ T> x x = 1Z X PI A x = span(Il(j ;i ,j ;i ,) e {_ lj o ) i}x{-i,o,i} \ a i x jx )( a i x jJ)> where l a oo) e {|0), |i)}, 
and the rest of \a ixjx ) G {|+), |-)}, for (i x ,j x ) ^ (0,0). For all y G N \ (0,0), V x „ y>x = 
lZ x _ y nA x = C <S)fce{-i o i}x{-i o 1} ^- Then it is clear that dim spa.n(Jl y£ ^V x _ yjX ) = 2 9 , while 
dim A x = 4 9 . Theorem 13. 161 then implies that this QCA is also not a QLGA. 

4.3. Examples of QCA that show effect of quantization on neighborhood. Consider 
a one- dimensional classical CA as in Figure [6] in which a cell is 2 bits. We consider the scheme 
in which a cell x is given by a pair (a x , b x ). The neighborhood is Af = {0, 1}, and the update 
rule /i: 

jj : a x i—)- a x 

b x ^-b x + a x + a x+1 

Let us quantize this CA to get a QCA. A cell x of the QCA consists of 2 qubits: \a x ) <S> \b x ), 
where a x ,b x G Z/(2). The quiescent state is |go) = |0) <S> |0). 
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Define the global evolution R on "He through its action by a controlled-NOT operation on 
the basis of its neighborhood. 

R : \a x ) H> \a x ) 

\b x ) H> \b x + a x + a x+ i) 

The neighborhood looks like A/+ = {0, 1} as expected. 
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Figure 6. One-dimensional CA whose neighborhood changes after quantization 

But we can equivalently look at the qubits in state: \a x ) ® |(— l) bx ), where a x , b x G Z/(2). 
Since |+) and |— ) are characters for the shift representation, R can be given as: 

R: \a x ) ^ (-l) axib *- 1+b * ] \a x ) 

K-i) 6 *) ^ |(-i) 6 *) 

Now the neighborhood looks like A/"- = { — 1,0}. In fact, the neighborhood is M = A/"_UA/+ = 
{ — 1,0,1}. This is a consequence of the well-known fact that the control and target of a 
controlled-NOT operation in the basis |+), |— ) are reversed as compared with the logical 
basis. Hence, while controlled-NOT operations can only propagate classical information in 
one direction (from control to target), they may propagate quantum information in both 
directions. 

As before, we compute the subalgebras V x - VjX = !Z x - y D A x , where y G Af. V X)X = 
H X V\A X = span({|a :I .)(a :r |(g)|6 x )(6 x .|}), where \a x ) G {|0), |1)}, and \b x ) G {|+), |-)}. V X - 1)X = 
TZ x ^i n A x = CI ® I. V X+1)X = TZ x + i n Ax = CI <S> I. Since dim span(II l , e {_i ) o ) i}X» a ._j /(X ) = 4, 
while dim A x = 16, Theorem 13.161 implies that this QCA is again not a QLGA. 

The change in neighborhood on quantization for this cellular automaton should lead us to 
question the original definition. Is there some redefinition of the cells that can maintain the 
same neighborhood for classical and quantum automata, and also yield a QLGA? We can 
in fact take the classical CA of Figure [6] and change it slightly by allowing a shift of bits, as 
shown in Figure [7] (the rearrangement of the add operations does not affect the original CA). 
Then, after quantization, we obtain a QLGA from it. We accomplish that because the shift 
of bits corresponds to a shift of tensor factors (advection). This QLGA has a neighborhood 
AT = {0, 1}, i.e., is a radius-1/2 QLGA. 
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Figure 7. One-dimensional radius- 1/2 QLGA from permutation of bits in 
the CA of Figure i 



5. Conclusions 

We have given a local algebraic condition for a quantum cellular automaton to be a 
quantum lattice gas. This result classifies QLGA as a subset of QCA, and resolves the 
question of which QCA are equivalent to a circuit with a "brick structure". This question 
has been open since Arrighi et al. [3] produced an example showing that QLGA and QCA 
are distinct classes, contrary to the work of [45]. The question of classifying the remaining 
QCA that are not QLGA remains open. This open question motivates the reexamination 
of the definition of QCAs through local rules [13,33,55]. The techniques introduced in this 
paper allow us to obtain rigorous results for models defined on the Hilbert space of finite 
configurations. The definition of this space necessarily requires the definition of products 
with countably infinitely many terms [23,51], just as is required for the constructions for 
WQCA defined through local rules [13,33,55]. One might hope that the techniques that are 
used in the present paper may also have utility in deciding which quantum automata defined 
through local rules are also causal, and therefore true QCA. 
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Appendix A. Some background and results on representations of 

ASSOCIATIVE ALGEBRAS 

We give a brief primer on finite-dimensional representations of associative algebras relevant 
to this paper. The material in this section, except for Corollary IA.6I and Theorem IA.81 is 
drawn from the texts by Goodman and Wallach [21, 22] and condensed for our purposes. 
Vector spaces in this section are assumed to be finite-dimensional and complex, unless 
otherwise stated. A representation of an associative algebra A is a pair (7, V) such that 
7 : A — > End(V) is an associative algebra isomorphism. Since a representation (7, V) of A 
makes V an ^4-module under the action: av = j(a)v, for a G A, and v G V, we can call V an 
.4- module instead (with 7 implied by the context). If U is a linear subspace of V such that 
j(o)U C U for all a G A, then U is invariant under the representation. A representation 
(7, V) is irreducible if the only invariant subspaces are {0} and V. 

Let (7, V), ([/,, W) be two representations of an associative algebra A. Let Hom(V, W) 
be the space of C-linear maps from V to W. Denote by Hom^V,!^) the set of all T G 
Hom(V, W) such that T^y(a) = fi(a)T for all a G A. Such a map is called an intertwining 
operator between the two representations. Two representations (7, V), (/i, W) are equivalent 
if there exists an invertible intertwining operator between the two representations. 

Lemma A.l (Schur's Lemma, Lemma 4.1.4, pg. 180, Goodman and Wallach [22]). Let 

(7, V) and (/i, W) be irreducible representations of an associative algebra A. Assume that V 
and W have countable dimension over C. Then 

1, if{>y,V)^1ji,W) 
0, otherwise 



dim Hornby, W) 



Let (7, V) be a finite-dimensional representation of A. Let A be the set of all equivalence 
classes of finite-dimensional irreducible representations of A. For each A G A, fix a model 
(ti x ,V x ). Let: 

U x = Rom A (V x , V) (37) 

For each A G A, U x is called a multiplicity space. 
Define the map: 

S x : U x <g> V x — > V (38) 
u x ®v x ^ u x (v x ) 

Then S\ is an intertwining operator with U X ®V X an ^4-module under the action a.(u®v) = 
u ® (av) for a G A. 

For A G A, define the X-isotypic component: 

V® : J2 W 

WCV:W~\ 

A finite-dimensional ^4-module V is called completely reducible if for every invariant sub- 
space W C V there is a complementary invariant subspace U such that V = W © U. The 
following result (Proposition 4.1.11, pg. 183, in [22]) provides equivalent characterizations 
of completely reducible representations. 

Proposition A. 2. Lei (7, V) be a finite-dimensional representation of A. The following are 
equivalent: 
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(i) V is completely reducible. 

(ii) V = V\ © • • • © V s with each Vi an irreducible A-module. 
(Hi) V = W\ + ■ • • + Wd with each Wj an irreducible A-module. 

Let (7, V) be a finite-dimensional representation of A, such that V is a completely re- 
ducible ^4-module. Then it is a result from the representation theory of associative algebras 
(Goodman and Wallach [21], Proposition 3.1.6, pg. 120), that the map S\ in (}3"H|) is an 
.4-module isomorphism: 

U x © V x ^ V [x] (39) 

Further: 

V = 

XeA 

This is called the primary decomposition of V. Thus if we extend the maps S\ as follows: 

5 7 = 0S A : 0[/ A <g>y A — ► V (40) 
Ae.4 Ae^ 

XeA XeA 

then Sry in (140 p is an isomorphism of ^.-modules (where V is an ^.-module via 7). Then we 
have the following isomorphism of algebras. 

S^jia) S' 7 = 0I t/ A®7r A (a) 
Ae^ 

where a G A, U x is as in ( 157)) . and (7r A , K a ) is a model for the class A. 

We state the Double Commutant Theorem, in Goodman and Wallach [22] pg. 184, special- 
ized here for a finite-dimensional vector space. Let V be a finite-dimensional vector space. 
Define, for any subset S C End(V), the commutant of S: 

Comm(S) = {y G End(V) : ys = sy Vs G S} (41) 

Theorem A. 3 (Double Commutant, Theorem 4.1.13 [22]). Suppose A C End(V) is an 
subalgebra containing the identity operator. Let V be a completely reducible A-module. Set 
B = Comm(^4). Then A = Comm(i5). 

Denote by Spec (7) the set of irreducible representations that occur in the primary decom- 
position of V. For each A G Spec(7), the multiplicity space U x = Hohi^Va, V) in ( l4"U)) is a 
Comm(^4)-module under the action: 

b.u = bu (42) 

where b G Comm(^4), and u G U x , and bu : V x uV bV is the composition (left multiplication 
by b) map. 

We adapt some results from Goodman and Wallach [22], Section 4.2.1 (General Duality 
Theorem). Except for some substitutions, these results are taken exactly as in the said 
reference, including the proofs. We obtain from the originals, derived in the setting of group 
representations, the versions relevant to subalgebras of End(V). Thus we get the associative 
algebra counterparts of the Duality Theorem, Theorem 4.2.1 [22], and Corollary 4.2.4 in [22]. 
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Theorem A. 4 (Duality, Theorem 4.2.1 [22]). Suppose A C End(V) is an subalgebra con- 
taining the identity operator. Let V be a completely reducible A-module. Each multiplicity 
space U x in ( |4(jp is an irreducible Comm(A) -module. Furthermore, if A, fi G A, and U x = £/ M 
as Comm(A) -module, then A = fi. 

Proof. We first prove that the action of Comm(^4) on U x is irreducible. Let T G U x be 
nonzero. Given another nonzero element 5 G U x we need to find r G Comm(.A) such 
that rT = 5. Let X = TV X and K = SV X . Then by Schur's Lemma (Lemma lA.ip . X 
and Y are isomorphic ^.-modules of class A. Thus there exists u G Comm(^l) such that 
uT : V x — > SV X is an ^4-module isomorphism. Schur's Lemma implies that there exists 
c G C such that cuT = 5, so we may take r = cu. 

We now show that if A ^ \i then U x and f/ M are inequivalent modules for Comm(^4). 
Suppose 

: U x — >U>* 

is an intertwining operator for the action of Comm(.A). Let T G U x be nonzero and set 
5 = 4>{T). We want to show that 5 — 0. Set {7 = TK A + 5V\ Then since we are assuming 
A 7^ fi, the sum is direct. Let p : £/ — > SV^ be the corresponding projection. Then 
p G Comm(^4). Since pT = 0, we have: 

= cj)(pT) = p<j>(T) =pS = S, 

which proves that = 0. □ 

Corollary A. 5 (Corollary 4.2.4 [22]). Let B = Comm(^l). Then V is a completely reducible 
B-module. Furthermore, the following hold: 

(i) Suppose for every A G Spec(7) there is given an operator T\ G End(K A ). Then there 
exists T G A that acts byl®T\ on the X-summand in the decomposition (I40p . 

(ii) Let T G A fl B (the center of A). Then T is diagonalized by (140]) and acts by scalar 
T(A) G C on U x ®V X . Conversely, given any complex valued function f on Spec(7) ; 
there exists T G A fl B such that T(A) = /(A). 

Proof. Since V is the direct sum of i3-invariant irreducible subspaces by Theorem IA.41 it is 
a completely reducible £>-module by Proposition IA.21 We now prove the other assertions. 

(i) Let T G End(V) be the operator that acts by I <E> T\ on the A summand. Then T G 
Comm(i3), and hence T G A by the Double Commutant Theorem (Theorem IA.3H . 

(ii) Each summand in fT40|) is invariant under T, and the action of T on the A summand 
is by an operator of the form R\<S>I = I® 5a with R\ G End(U A ), and 5a G End(V A ). 
Such an operator must be a scalar multiple of the identity operator. The converse 
follows from (i) 

□ 

An algebra B is simple if the only two sided ideals in B are and B. An algebra A 
is a semisimple algebra if it is a finite direct sum of simple algebras. This is equivalent, 
by Wedderburn's Theorem (Goodman and Wallach [21], Theorem 3.2.1, pg. 128), to the 
statement that A is isomorphic to a finite sum of matrix algebras. 

For an algebra A C End(V^), such that V is a completely reducible ^4-module, let Spec(^4.) 
be the set of irreducible representations that occur in the primary decomposition of V. 
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Corollary A. 6. Let A C End(V^) be a subalgebra containing the identity operator, such 
that V is a completely reducible A-module under the action of End(V) restricted to A. Let 
B = Comm(^4). Then there is an A-module isomorphism: 

S V :($U X ®V X —>V (43) 
XeA 

XeA XeA 

where V x is an irreducible module for the class X G A, and U x = }iom^(V x , V). Under this 
isomorphism: 

A = Q)lux(g)End(V x ) 
XeA 

and: 

B = 0End([/ A ) ®I v a 

where the action of B on each U x is given by (|42|) . 
In particular, A is semisimple. 

Proof. V is a completely reducible ^.-module. This implies there is an ^4-module isomorphism 
given by (140]) : 

S v : U x ® V x — > V 

XeA 

^2 uX ® yX ^ 5^ uX ( vX ) 

XeA XeA 

where V x is an irreducible module for the class XeA, and U x = Hom^(V A , V). Then by def- 
inition of A action (as a subalgerba of End(V^)) on V, under this isomorphism, Corollary IA. 51 
implies: 

.4= V®End(y A ) 

Specpl) 

Since V is finite-dimensional, the above implies that Spec (A) is a finite set. This in turn 
implies that A is semisimple as it is isomorphic to a finite direct sum of finite dimensional 
matrix algebras. Then a standard theorem on semisimple algebras, Proposition 3.3.1 in [21], 
implies that each irreducible representations of A is equivalent to some element of Spec(^4). 
Hence A = Spec(^4). Therefore: 

A = V g> End{V x ) 
XeA 

By Duality Theorem (Theorem IA.4j) . under Sv'- 

B 0End([/ A ) ®I v x 
XeA 

where the action of B on each U x is given by (j42p . □ 
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We state the semisimple algebra version of the Double Commutant Theorem, in Goodman 
and Wallach [21] pg. 137. 

Theorem A. 7 (Double Commutant Theorem for Semisimple Algebras). Suppose A C 
End(V) is a semisimple subalgebra containing the identity operator. Then the algebra B = 
Comm(.4.) is semisimple and A = Comm(i3). Furthermore, there exists an A-module iso- 
morphism: 

r 

S A : Uj ® Vj — > V 

r r 

j=i i=i 

where Vj is an irreducible A-module, and Uj = Hom_4(V 7 -, V). Under this isomorphism: 

r 

^S0I p .<g)End(^) 
i=i 

and: 

r 

B^0End(C/ i )«)Iy. 
i=i 

Next we consider a general Hilbert space, and state the remarkable von Neumann Density 
Theorem. The proof of the von Neumann Density Theorem can be found in [26,53]. Consider 
a Hilbert space % and the space, B(l-L), of bounded linear operators of %. Define, for any 
subset S C B(%), the commutant of S: 

Comm(5) = {y G B(H) : ys = sy Vs G S} 

Theorem A. 8 (von Neumann Density Theorem). Let A C B(H) be a self-adjoint algebra of 
operators containing the identity operator, then the weak and strong closures of A in B(7i) 
is Comm(Comm(^4)). 



Appendix B. Some results on representations of self-adjoint algebras 

We state a few results on representations of self-adjoint algebras. Let V be a finite- 
dimensional Hilbert space, with an inner product (•[•), and A C End(V^) a subalgebra. Let 
the adjoint of an algebra A C End(K) be AJ = {A^ : A G A}. An algebra A is self-adjoint 
if A = A*. From the results in this section, Lemma [B. 11 and Proposition IB. 21 appear in [42] 
(pg. 145). 

Lemma B.l. Let A C End(V) be a self-adjoint subalgebra. If W G V is an A-invariant 
subspace, then W L = {v G V : (v\w) = OVw G W} is A^ -invariant. 

Proof. Let w G W, v G W 1 , A G A. Then A ] G A, which implies A^w G W =^ 
(Av\w) = (v\A^w) = =^ Ah e W ± . □ 



Proposition B.2. Let A C End(V) be a self-adjoint subalgebra. Then V is an orthogonal 
direct sum of irreducible A-modules. In particular, V is a completely reducible A-module. 
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Proof. Let W C V be an ^.-invariant subspace of minimal dimension. Then it is by definition 
irreducible. Since A = A' , by Lemma IB.lt V = W @ W 1 - is an orthogonal direct sum of 
^4-modules. The conclusion follows by induction on dimension. □ 

Corollary B.3. Let A C End(V) be a self-adjoint subalgebra containing the identity opera- 
tor. Then A is semisimple. 

Proof. By Proposition IB.21 V is a completely reducible ^4-module. Then Corollary IA.6I 
implies that A is semisimple. □ 

Corollary B.4. Let A C End(V) be a self-adjoint subalgebra. Then there exist for every 
A G A, a set of intertwining operators {u x } C U x , such that the X-isotypic component can 
be written as an orthogonal direct sum: = (&jU X (V x ). 

Proof. By Proposition IB.21 V is a completely reducible A- module. Thus by (1391) . S\ in (1381) 
is an ^4-module isomorphism: 

S x : U x g> V x — ► V [x] 
Again by Proposition \B.2\ the A-isotypic component, 



vn = w s 

WjCV:Wj~\ 

is an orthogonal direct sum. Choose Uj such that Uj(V x ) = Wj. □ 

Appendix C. Structural Reversibility 

We prove the Structural Reversibility (due to Arrighi et al [3]), in our context, to make 
the present paper self-contained. 

Theorem C.l (Structural Reversibility ). Let M : He — > "He be a unitary operator and 
M a neighborhood. Let V = {— k\k G A/"}. Then the following are equivalent: 

(i) M is causal relative to the neighborhood N '. 

(ii) For every operator A x local upon cell x, M^A X M is local upon J\f x . 
(Hi) Aft is causal relative to the neighborhood V. 

(iv) For every operator A x local upon cell x, MA X M^ is local upon V x . 

Before proving the theorem, we prove the following lemma: 

Lemma C.2. Denote the set of density operators on % c by fl. For A G B(7-L c ), define the 
map: 

f A -.n^c 

f A (p) = tr(Ap) 

Let DgZ" be a finite subset. Then: 

(i) A is local upon D if and only if for every pair p, p' G £1 satisfying p\ D = p\ , Ja{p) — 

Up')- 

(ii) p, p' G Q satisfy p\ D = p'*, if and only if for every A local upon D, /_a(p) = /^(p')- 
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Proof. For elements of B^Hc) that are local on some finite subset K C Z n , by abuse of 
notation, we call both B G Z {Z as in (FIU]) ) and l]^(B) (l k as in by £>. 

Suppose A is local upon D and p, p' G Q are such that p\ D = p'\. Then /U(p) = 



tr(Ap) = tr(Ap,J = tr(V| J = MV) = /a(p')- 

For the converse, suppose A is not local upon D. We can choose an orthonormal basis 
{|%)} for <S) keD W. Then we can write (special case of ([6])): 

A = ^2 \ Vl )( Vm \ ® Ai,m 
(l,m) 

where {(I, m)} C {1, . . . , dw} x {1, • • • , dw}, and A; iTn G B(7ic—) are non-zero. Since A is 
not local upon D, not all the A^ m are multiples of the identity, ILj, on There are two 

cases to consider: 

(1) There is an Aij for some I such that Aij is not a multiple of %. Then there exist 
two unit vectors \x), \y) G Hc^, such that (x\Aij\x) ^ {y\A^i\y). Let p = |t>z)(t>;| ® 
|x)(x|,p' = \vi)(vi\ g> |j/)(y|. Then p b = p' b whereas / A (p) 7^ /a(p')- 

(2) Every is a multiple of I-^. Then choose some A^ m which is not. Then there 
exist two unit vectors \x), \y) G such that (x\Ai tm \x) ^ (y\Ai ym \y). Let p = 
{\vi)(vi\ + \vt)(v m \)®\x)(x\,p' = {\vi)(vi\ + \vi)(v m \)®\y){y\. Then p\ D = p\ D whereas 
/a(p) + /a(p')- 

Suppose A is local upon D and p, p' G ft are such that p| D = p'| D . Then /a(p) = 



tr(Ap) = tr(Ap lD ) = tx{Ap\ D ) = tr(Ap') = f A (p>). 

For the converse, let p, p' G ft, such that for all A local upon Z), /U(p) = /a(p')- Consider 
A = |t>i)(i> TO | ® %, a basis element of operators local upon D. Then /U(p) = (v m \p\ D \vi) , 
/a(p') = (^mlp'ioh)- This implies by /a(p) = /a(p')> tna t (f m |p| D h) = (^m|p'| D k)- Since 
this is true for all elements of {\vi)(v m \ : 1 < l,m < dw}, i-e., for a set of basis elements of 
operators local upon D, then by linearity, p\ D = p'<. 

□ 

Proof of Theorem I C. 1\ 
Note: 

tr(MpM f A) = tr(pM f AM) : Vp G ft (44) 



(i) =r- (ii) Let M be causal. Suppose p and p' is a pair such that pl^. = p'\j\f x , and 
A is local upon x. Causality implies tr(MpM^A) = tr(Mp'M t A). By dSJ), ti(pM^AM) = 



tr(p'MUM), so by Lemma E^fM AM is local upon M x . 



ii. 



Suppose A is local upon x, p and p' be a pair of density operators such that 
p\^~^ p'\ N ~uc\ MUM is local upon Af x . Then tr(pMUM) = tr(p'ML4M). By (gU), 



tr(MpMU) = tr(Mp'MU), so by Lemma Op)] MpM\. = Mp'M^ 



IV 



Suppose A x is local upon x. Let A y be an operator local upon y ^ V x . 
ThenM t A y MTTlocal upon My. But x ^ A/^, so M^AyM commutes with A x . But then 
A y = MM^ AyMM^ commutes with MA X M\ by the definition of conjugation. This applies 
to every A y G A y , hence MA X M^ commutes with A y . We can express MA X M^ as a finite 
sum (as in ([6])): 



MA X M^ = ^e l ®E l 
i 



(45) 
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where {e{\ C End((£) keVx W) ia a linearly independent set, and E\ G B(7i C —) (the space of 
bounded linear operators on the co-V x space, He , defined in Definition 13. 1[ with D = V x ). 

V x 

For any B G End((3) fcgVi W), B eg) Ip- (where Ipj is the identity on "He—) also commutes 
with A y , This implies (B ® Iy-){MA X M^) commutes with A y . Using the standard commu- 
tator notation [., .] ([X,Y]=XY-YX), we write: 

[A y , (B ® I W )(MA X M^)} = J2( Be i) ® \K Ei\ = 



We can trace out the finite dimensional part. Let the operator thus obtained on He be 

Vx 

KB): 

X(B) = J^tr (Be^Ay^t] = ^tr (Bei)[A y ,Ei] = 
i i 

Since trace is a degenerate bilinear form on End(0 feeV ^ W), and {e;} is a linearly independent 
set, we can choose a set {B{\ (indexed by the same set {1} that indexes {e;}) such that 
tr (B k ei) = S k) i- 

\{B k ) = [Ay, E k ) = (46) 
The subalgebra of BCHq ) generated by A v , y 4. V x acts irreducibly on the space He 

Vx _ V x 

(by the argument in the proof of Theorem 13.7} i.e., it acts transitively on a dense sub- 
space, span(Cp^), of He—)- By f T4"6"j) . E k commutes with this subalgebra. Schur's Lemma 
(Lemma lA.ip then implies that E k = ?7fcly- for some 7] k G C. This is true for every k G {/}. 
BygSJ: 

MA X M^ = (J2 Viei) ® % = a Vx ® % 



for some ay x G End((3) fceV ^0- Hence MA X M^ is local upon V x . 
Same argument as for the proof of 



iv , 



n. 



the roles of Af x and V x , M and ML 

Same arguments as for the proof of 



iv 



except interchanging 



in. 



iv 



with V, M x with V x , and M with ML 



ii. 



except replacing M 
□ 



Appendix D. Proof of Lemma 13.151 requiring quiescent symbol \q ) to be an 

INVARIANT OF F 



Lemma (Lemma I3.15p . F is defined on C as a map F : C — > T-Lg, and can be extended to 
H c as unitary and translation-invariant operator, if and only if F has {%) as an invariant 
(an eigenvector with eigenvalue one): F\q ) = \q ). 



Proof of Lemma VJ . 1 5\ Let us assume that F is defined on C as a map F : C — > H^, and 
can be extended to H^ as a unitary and translation-invariant operator. Then Lemma 13.51 
implies: 



F( <g) |gb» 



F\q ) 



\Qo) 



J't- 



for some 0o G R. This implies, F\q ) = \q ). 

For the converse, assume F\q ) = \q ). Then for every element |c) = ® x6Z n \c x ) G C, 
under the formal definition of F in (1261) . F\c) = ® xeZ n F\c x ) G Tig. Furthermore, as 
a map F : C — > %£, F preserves the inner product (since F is a unitary operator on 
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W = <S>zeAf an< ^ * s translation- invariant, by the same definition. Then F can be extended 
to a unitary and translation-invariant operator on T-Lq. □ 



